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We make a detailed calculation of the amplitude of isocurvature perturbations arising from inflation- 
ary models in which the cold dark matter is represented by a scalar field which acquires perturbations 
during inflation. We use this to compute the normalization to large-angle microwave background 
anisotropies. Unlike the case of adiabatic perturbations, the normalization to COBE fixes the spec- 
tral index of the perturbations; if adiabatic perturbations are negligible then rii so ~ 0.4. Such blue 
spectra are also favoured by other observational data. Although the pure isocurvature models are 
unlikely to adequately fit the entire observational data set, these results also have implications for 
models with mixed adiabatic and isocurvature perturbations. 
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I. INTRODUCTION 



There has been recent interest pHq] in inflationary 
models which give rise to isocurvature perturbations in 
the Universe, rather than the traditional adiabatic ones. 
The simplest way to achieve this is to assume that the 
material presently making up the cold dark matter takes 
the form of a massive scalar field, which during inflation 
received perturbations by the usual mechanism. With 
suitable subsequent evolution, these perturbations can 
be the dominant ones in the present Universe. 

The familiar inflationary calculation for adiabatic per- 
turbations is extremely simple (see Ref. jjj for a review), 
because once the perturbations have been stretched to 
wavelengths greater than the Hubble length they become 
constant, and no calculation is required until the much 
later epoch when they cross back inside the Hubble ra- 
dius, which for scales of interest is well after nucleosyn- 
thesis. By contrast, isocurvature perturbations (| [n]] 
are free to evolve on superhorizon scales, and the am- 
plitude at the present day depends on the details of the 
entire cosmological evolution from the time that they are 
formed. On the other hand, because all super-Hubble ra- 
dius perturbations evolve in the same way, the shape of 
the isocurvature perturbation spectrum is preserved dur- 
ing this evolution and so need not be calculated. In the 
literature, therefore, comparison with observation |3]]4|] 
has largely focussed on the spectral index of the pertur- 
bations, and it has been assumed that the model param- 
eters could always be tweaked to ensure the amplitude is 
correct. However, as we will see the constraints on ampli- 
tude and slope are intimately related in these scenarios. 
In this paper, we compute the Sachs- Wolfe contribution 
to large-angle microwave background anisotropies for a 
particular class of isocurvature models which give non- 
gaussian perturbations. Unlike the case of adiabatic per- 
turbations, the Sachs- Wolfe effect fixes the slope of the 
isocurvature perturbations, as well as their amplitude, 
because the slope determines the matter density. 



II. MODEL DESCRIPTION 

The class of models we discuss is characterized by 
the cold dark matter being a scalar field with a time- 
dependent mass, which apart from its quantum fluctua- 
tions sits in the vacuum state during the entire history of 
the Universe. Particular implementations of that general 
statement have been given by Linde and Mukhanov 
and by Peebles M. It is well known that a sufficiently 
massive scalar field 4>, with m^$> H where H is the Hub- 
ble parameter, oscillates with decaying amplitude such 
that its energy density p$ oc 1 /a 3 where a is the scale fac- 
tor, and further that perturbations in this field obey the 
same equations as perturbations in a pressureless fluid. 
Such a scalar field is therefore a viable cold dark mat- 
ter candidate, the classic example being the axion (see 
Ref. [[l0| for an early treatment and Ref. jnj for a re- 
view) . 

The mass of the scalar field in such a scenario must 
have quite a complicated evolution in order to generate 
a viable model. 



A. During inflation 

If the field is effectively massless, m <C H, its per- 
turbations will be scale-invariant. Such a spectrum is 
excluded by comparison of COBE data with galaxy clus- 
tering; what is needed is a spectrum with a tilt favour- 
ing short-scale perturbations (usually called a blue spec- 
trum). This can successfully be generated if m ~ H at 
the relevant stages of inflation jl[ ; the mass causes per- 
turbations to die away outside the horizon, and those 
which have been outside the longest receive the greatest 
suppression. However if m H during inflation, then 
perturbations will be highly suppressed with an nj SO = 3 
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spectrum,|j which is also unsatisfactory. One must there- 
fore tune the mass to be close to H, perhaps best done 
by generating an effective mass m 2 ff = aH 2 where a is 
some coupling hoped to be of order one. 

B. After inflation 

We are seeking to create isocurvature perturbations, of 
the type where the CDM perturbation is cancelled by the 
collective perturbation of all other species. In fact, the 
perturbations we are discussing are actually isothermal, 
since the fluctuations in cj> are uncorrelated with those 
in the inflaton field which is what ultimately determines 
the radiation density after reheating. These isothermal 
perturbations can be broken into an adiabatic and an 
isocurvature part. However, as long as the energy den- 
sity of the CDM field is totally subdominant to that of 
the inflaton at the end of inflation, the isocurvature per- 
turbation will be by far the dominant part M . 

Another important point is that the inflaton must 
not decay into the CDM field, which would destroy the 
isocurvature perturbations, and so must decay only into 
conventional particles. Therefore at the end of inflation 
the CDM field has much less energy density than every- 
thing else, whereas by the present it has to dominate the 
baryonic component (e.g. the present Universe may have 
^baryon — 0.04, ~ 0.3 and a cosmological constant). 
While the field has a significant mass then p^ cx 1/a 3 , as 
compared to the thermalized fluid which at early times 
goes as /9 ra d oc 1/a 4 . However it turns out that a more 
dramatic gain is required, which is achieved by having an 
epoch where the mass is small enough to be negligible. 
When this happens, the scalar field feels no force and is 
effectively fixed, leading to a constant enabling it to 
'catch up' with the conventional matter. 

C. At late times 

In order to preserve the successes of the standard cos- 
mology, the field preferably should already be acting as 
cold dark matter by the time of nucleosynthesis (though 
there may be no observational conflicts even if this be- 
haviour is delayed until sometime later). This requires 
that the mass is greater than the Hubble parameter at 
that time, and this guarantees that it will act this way 
when any astrophysically-interesting perturbations enter 
the horizon. Depending on whether one believes in a 
cosmological constant or not, the desired dark matter 



*We use conventions where ni so = corresponds to scale- 
invariance. Unfortunately the literature features at least three 
different conventions, where scale-invariance can correspond 
to an index of —3, or 1. 



density is between about 5 and 20 times the baryon den- 
sity, so the field must start to scale as 1/a 3 just when it 
reaches that relative density. This requires fine-tuning, 
but that's true of all dark matter models as they are all 
trying to explain the same coincidence of dark matter 
and baryon densities. The simplest scenario is that the 
mass is the bare mass of the field; the effective mass was 
large during inflation due to some coupling, and then the 
field evolved as effectively massless until H fell to be of 
order this bare mass. 

An interesting feature of this model is that the pertur- 
bations are crucial to define the background. In normal 
models there is a homogeneous background density /Ocdm 
about which the perturbations are defined, but here the 
homogeneous background is <p — =>■ p^ = 0. Instead, 
the mean density of the dark matter is generated by the 
perturbations, taking advantage of the density being pro- 
portional to the square of the field perturbation. At each 
point in space the field <fi is oscillating with some am- 
plitude. If we define ^(x) to be the maximum of the 
oscillation at point x, then the energy density is 

W (x) = lm 2 <D 2 (x). (1) 

and the mean density is 

p = im 2 (<l> 2 (x))=m 2 (^(x)) 

= m 2 f^Vm^ • (2) 

In that last expression, V is the power spectrum of 0, 
defined as 

= |^0k) , (3) 

where </>k are the Fourier components of </>(x) and k = |k|. 
We have included an upper cut-off A: max , intended to be 
the shortest scale on which there are perturbations and 
estimated as the Hubble scale at the end of inflation. 

The background dark matter density therefore has 
to be obtained self-consistently out of the perturbation 
calculation. Were the spectrum flat, this cannot be 
achieved; we need the power spectrum to rise sharply 
to short scales, as then the correct density can be ob- 
tained by suitable tuning of fc max and m to give the ob- 
served value. Calculationally this is also desirable, as in 
considering the large-scale perturbations one can imag- 
ine smoothing out the small-scale perturbations to ob- 
tain a genuinely smooth background. It is interesting, 
though certainly not compelling, that the blue spectrum 
is favoured both in order to obtain the correct dark mat- 
ter density and to match large-scale structure observa- 
tions. 

One should also note that the statistics of the per- 
turbations are neither gaussian nor strictly speaking chi- 
squared, though they are related to the latter. It is the 
total density p^ which is chi-squared distributed, not the 
density contrast <5(x) = [p^(x) — p^/p^. 
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III. BACKGROUND AND PERTURBATION 
EVOLUTION 

The complete model needed to evaluate the amplitude 
of perturbations therefore has two inputs. The first is the 
evolution of the Hubble parameter H(t). During infla- 
tion this is governed by the inflaton field, which is some 
other scalar field whose energy is dominating over the 
CDM field. After inflation, it is determined by whatever 
the matter content of the Universe is, and we shall be 
assuming the standard, simplest evolution of radiation 
domination after inflation giving way to matter domina- 
tion. One could include a cosmological constant too but 
we will ignore that minor complication. The second in- 
put is the time-dependent mass of the CDM field. As 
discussed, this will be of order H during the late stages 
of inflation, then at the end of inflation drops to a small 
constant value which is initially negligible. Then even- 
tually the Hubble parameter drops below this value and 
the field begins truly to behave as CDM. 

The perturbations need to be tracked at various stages. 
First the amplitude arising from quantum fluctuations 
needs to be computed. Subsequently, the perturbation 
is outside the Hubble radius, but still decays as it has a 
significant mass m ~ H . At the end of inflation the mass 
becomes negligible and the field freezes out, before even- 
tually the mass becomes important and the field begins 
to evolve again. For all the earlier stages, it is simplest 
to define the evolution in terms of the scalar field itself. 
However, at the last stage we have exactly a standard 
isocurvature CDM model, and standard results (for ex- 
ample for the Sachs-Wolfe effect) can be applied. These 
are usually expressed not in terms of the field, but rather 
in terms of the entropy perturbation S defined later on. 

To begin with we need to study perturbations in the 
CDM field which is responsible for generating isocurva- 
ture fluctuations. The background evolution is governed 
by the inflaton field, which dictates that the Hubble ex- 
pansion during inflation is H 2 ~ V/3Mpj, where V is 
the potential for the inflaton field and Mpi = v/1/ 8ttG 
is the reduced Planck mass. As explained earlier, we as- 
sume that energy density of the CDM field <fi is subdom- 
inant compared to that of the inflaton during the slow- 
roll phase. We are interested in studying the fluctuations 
generated in the cj> field. The fluctuations in the inflaton 
field cause perturbations in the curvature giving rise to 
adiabatic perturbations, but such perturbations are to be 
arranged to be negligible compared to the isocurvature 
ones. 

There are various sources which would produce an ef- 
fective mass to the CDM field 0. One possibility is the 
existence of a time-varying coupling of the CDM to the 
inflaton field. During inflation such a coupling can be ad- 
justed to give an effective mass to the CDM which is pro- 
portional to the Hubble parameter. Another possibility 
is the existence of a non-minimal coupling of the CDM to 
gravity. Such a term appears as a correction to the mass 



of the CDM and does not alter the dynamics of the over- 
all evolution as the CDM field is oscillating at the bottom 
of the potential = 0. Such a coupling gives a correction 
to the mass of the CDM field: m 2 (H) = m 2 + f 2^H 2 @. 
For our purposes, it does not matter what the source of 
the correction is as long as it fixes the effective mass to 
be of the order of the Hubble parameter. During infla- 
tion the presence of such a term provides the tilt in the 
power spectrum and just after the end of inflation it must 
switch off to avoid the CDM field dominating too early. 

To study the perturbations we need the perturbation 
equation for the CDM 



3H6cf> k 



«2 



effW 



k =0, 



(4) 



where the effective mass m c g(t) is 



m 2 eS (t) =m 2 +a(t)H 2 , (5) 

where m is the bare mass term of the CDM field. We 
assume that m 2 <C a(t)H 2 during inflation. The solution 
of Eq. (B) for k < aH is well known [O] 
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If the ratio of effective mass m c g(t) and H is constant 
during the entire period of inflation, and assuming H is 
suitably slowly varying, then at the end of inflation the 
solution for the power spectrum for the modes which are 
already outside the horizon is pj 



2tt- 



2tt 2 



H 



Hm c 



2a/3 



(7) 



where the subscript 'e' denotes the end of inflation. 

This final expression, Eq. (Q), is strictly correct only 
at the end of inflation, and we have replaced m c ff with 
aH 2 , where a is a constant less than one. An alter- 
native physical scenario would be if m e ff were constant 
rather than a, in which case the tilt will pick up an extra 
scale-dependence, depending on the details of the explicit 
model. We will only consider the situation when a is a 
constant. Even in that case, an extra scale-dependence 
of the tilt arises because of variation of H during infla- 
tion, but this contribution will be insignificant compared 
to the tilt induced by the mass. 

After inflation ends, the effective mass must drop pre- 
cipitously from its value during inflation to the bare mass 
m which is then much less than H . Initially then the 
mass can be neglected, and Eq. (|4|) has a simple solution 
on large scales: Stfik = constant. As a result the power 
spectrum, Eq. (0), remains unchanged. 

This simple situation does not hold for ever, as the 
Hubble parameter drops as a~ 2 during radiation era, 
where a is the scale factor of the Universe. Eventually the 
Hubble parameter becomes equal to the bare mass of <p 
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and henceforth the mass of the (j> field cannot be neglected 
in Eq. (0). From this time onwards, the CDM field acts 
like a true massive CDM field, and the amplitude of the 
perturbations falls as a~ 3 / 2 during all subsequent phases 
of the Universe. Noting that H(t) = H c a 2 /a 2 during ra- 
diation domination, we find that at any time after mass 
turn on the spectrum of the fluctuations is given by 



2a I Z 



(8) 



Preferably this transition should take place before the 
synthesis of light elements, when the ratio between the 
photon and the baryon is already a fixed constant, though 
compatibility with observations may still be maintained 
even if it happens some time afterwards. This expression 
for the spectrum continues to hold at all later epochs, 
including matter domination. 

After mass turn-on, it is useful to consider an alterna- 
tive quantity, the entropy perturbation S given by 



S 



3 Sp r 

4 p r 



(9) 



where the subscript V denotes radiation. Because the 
perturbations are actually isothermal (albeit in a limit 
where isothermal and isocurvature are almost the same 
thing), we have 5p r = and this can be written as 



5fj>_ _ P4>(x) - p<t> _ S(t) 2 (x) - (84> 2 {x)) 



(10) 



where p^ is the mean energy density and it is defined as 
the spatial average of the field fluctuations. p<f>{x) de- 
fines the total energy density of the field (f> at any given 
time and it is determined using Eq. (§) from the field 
fluctuations themselves. The advantage of using the en- 
tropy perturbation is that it is a standard result p2|j7[ 
that on super-horizon scales it continues to be conserved 
even when the radiation is no longer dominant (at which 
point the radiation will acquire perturbations from the 
influence of the CDM field fluctuations on the expansion 
rate), and this conservation law can be used to evolve the 
perturbations up to decoupling to compute the Sachs- 
Wolfe effect. 

Now we have all the necessary formulae to determine 
the Sachs-Wolfe effect for the isocurvature fluctuations 
of the 6 field. 



IV. COMPARISON WITH OBSERVATIONS 

A. Normalization from the Sachs— Wolfe effect 

The anisotropy due to the entropy perturbation is 
jiven by the well-known formula 0,|T3|,|7| 



AT(x) 
T 



2S_ 
~5~ 



(11) 



Followingthe same steps as for the adiabatic case (see 
e.g. Rcf. 0), the spectrum of the cmb, denoted by Q, is 
given by 



C.= 



16tt 
~25~ 



dk 
~~k 



2k 
H a a 



(12) 



where V s {k)S(k-k') = (S£Sv)k 3 /2tt 2 , and j e is a spher- 
ical Bessel function. We have taken the distance to the 
last-scattering surface as 2/ciqHo, where the subscript 
corresponds to the present value. This formula is strictly 
true only for a critical-density universe, but could be 
readily generalized to low-density. 

From Eq. (|l^) , except for the k = mode we have 
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Hence 



Ps(k) 



(k + 0) . 



(13) 



(14) 



To evaluate Ci we need to compute the power spectrum 
of Sip 2 . Its calculation is quite subtle, and we carry it out 
in the Appendix. The result, for the range of tilts we are 
interested in, is 



I(2a/3) /a e \ 6 / H, 



Air 4 



1 - 2a/3 \H e a 



4a/3 



(15) 



where I(2a/3) is a constant described in the Appendix, 
which is of order 10 for the values of a we will be inter- 
ested in. Here H e is the Hubble parameter at the end of 
inflation. 

Substituting Eq. ( |l5| ) into Eq. jl2]), we get the follow- 
ing expression 



C t = 



4/(2a/3)m 4 ( a 
1 

(ifeae) 4 "/ 3 Jo 



4 ' - \ - 6 'H n \ 3 H 4 



1 - 2a/3 



Eq. (^) can be further simplified by performing the in- 
tegration 



t(! + i)c t 



2 1(2a/3) m 4 H 4 a 6 c ff c 3 / H a„ 



2^p 2 (1 - 2a/3) 



\2H C 



4a/3 



^ r[(2-4a/3)/2]r[^+2a/3] ) 

2 [ J r[(3-4a/3)/2]r[^+(2-2a/3)]J ' 1 ' 



In fact, as long as £ 3> 1 and £ 3> 2a/3, the complicated 
term within the braces can be approximated simply by 

£^/3 ( gee e g Ref 

Although in this expression the quantities p^ and a 
were given at the time of decoupling, since p^ oc 1/a 3 
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they can in fact be evaluated at any epoch, including the 
present. 

Next we can estimate by using Eqs. (||) and (||), 
obtaining 



3m 2 
4n 2 a 



3/2 



H, 



2a/3 



(18) 



Substituting in Eq. (|l7|), and taking fc max ~ H c a c 
since we assume that the shortest scale perturbations 
generated are at the horizon scale at the end of infla- 
tion, one gets a simple relation for a multipolc of order 
£, namely 



e{e + i)c t ~— l _ 2a/z 



2iJ n et n 



4a/3 



(19) 



The COBE observations juj give the 10th multipole as 
1)Q - 7 x 10~ 10 . Taking the ratio of Hubble scales 
at the present and at the end of inflation to be given by 
the usual estimate In {H c a c / Hqclo) w 50, this gives 



4a 



0.43, 



(20) 



where nj SO refers to the slope of Vs- 

Predominantly, then, the Sachs- Wolfe effect for isocur- 
vature perturbations fixes the slope of the perturbations, 
rather than the amplitude as one finds for adiabatic 
perturbations. The reason is that the slope, favouring 
short wavelengths, is what determines the matter den- 
sity, and for the observed perturbations S ~ 10~ 5 to 
be achieved, the tilt must be steep enough to favour 
short-scale perturbations in determining the matter den- 
sity against which the perturbations are observed. Were 
the spectrum flatter, then the background density would 
not be high enough, whereas if it were steeper than the 
large-scale perturbations would not give enough density 
contrast against the mean background derived from the 
short scales. This argument was given in Ref. §; our 
calculation had given a precise quantification of the con- 
straint. 

It is intriguing that the slope required to give the cor- 
rect magnitude of the Sachs-Wolfe effect is also in bet- 
ter agreement as to the slope of the large-angle Ci than 
a flat spectrum, although it seems unlikely that a pure 
isocurvature model will be able to fit the complete set of 
observational data flla,01 . 



B. The CDM matter density 

We can constrain the bare mass to of the CDM from 
the present density 



Pcrit — 3Hq Mp! , 



(21) 



where we continue to assume that we live in a critical- 
density universe. With the help of Eq. ( |l8| ) we can derive 
a simple expression for the CDM mass 



4 H o 



(22) 



where we took /c max = H c a c .To simplify this, we model 
the Universe as matter-dominated during reheating until 
Hrh, followed by the usual cosmology. During matter 
domination a oc H~ 2 / 3 and during radiation domination 
a oc H^ 1 / 2 , leading to 



to = l^WM*^ 
HI 



<^eq ^rh 



167r 4 a 2 M^ 



H 



cq 



ffrh Hf 



(23) 



For critical density the Hubble parameter at matter- 
radiation equality is H eq = 5.25 x 10 6 h 3 Hg where Hq = 
2.1/i x 10" 42 GeV (with 0.5 < h < 0.75) is the present 
Hubble parameter. Putting in characteristic values leads 
to 



lGeV 



44a g, 



-1/2 , 4 



/l0 13 GeV\ /l0 10 GeV 



(24) 



where g* ~ 100 is the number of particle species at re- 
heating. 

The model appears compatible with observations for 
a reasonable range of parameters. Firstly, gravitational 
wave production requires H e < 10 14 GeV, or gravita- 
tional waves would give the dominant contribution to the 
COBE anisotropies (see e.g. @). Indeed, it will need to 
be some way less than this to make sure that adiabatic 
perturbations do not exceed the isocurvature ones which 
we presume to dominate. The mass to needs to be less 
than H e for the scenario to be consistent, and should not 
turn on too late if the standard cosmology is to be re- 
covered. Everything would be completely safe if turn on 
precedes nucleosynthesis, but may continue to be permit- 
ted even if turn on is later, since the dark matter den- 
sity would be negligible during nucleosynthesis. However 
since H nuc \ ~ 10 -24 GeV this is unlikely to be the situ- 
ation. Finally, if ones assumes supersymmetry then the 
reheat temperature is constrained by gravitino produc- 
tion from two-body scattering processes in the thermal 
bath. If the reheat temperature is too high then grav- 
itinos are produced with an abundance which destroys 
standard nucleosynthesis. To avoid this, it should satisfy 
be T rh < 10 10 GeV @. This implies reheating should 
be inefficient, though the problem is no different here 
to usual inflation models, and may indeed be less severe 
than some standard models as we can accommodate a 
low energy density at the end of inflation. 

Parameters exist which satisfy all these constraints. 



For example, H c ~ 10 12 GeV, T rh 
10 4 GeV is a suitable combination. 



10 9 GeV and 
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V. DISCUSSION AND CONCLUSIONS 

We have carried out a full calculation of the ampli- 
tude of microwave anisotropies within a class of isocur- 
vature inflation models which includes models by Linde 
and Mukhanov p| and be Peebles Q . We have not dwelt 
on the model-building issues of obtaining the appropriate 
time-dependence of the CDM field's mass, but under the 
assumption that it can be made to behave as required 
we have evaluated the model constraints. The key result 
is that the Sachs-Wolfe effect in these models primar- 
ily constrains the slope, rather than the amplitude, of 
the perturbations. Typically observational comparisons 
have assumed that the slope is a free parameter (e.g. 
without taking this additional restriction into account. 

A broad range of parameters appears capable of si- 
multaneously matching the COBE observations and the 
required CDM matter density. It is intriguing that the 
slope rtj so ~ 0.4 required to match the Sachs- Wolfe am- 
plitude is also in the right vicinity to have a chance of 
matching the slope of the large-angle anisotropy spec- 
trum. Nevertheless, it is probably true that such pure 
isocurvature scenarios are unable to fit the complete data 
set |l5|,|]]. There does however remain scope for mixed 
adiabatic and isocurvature scenarios |lf]|||l7],f|] , and our 
results can be applied, with minor modification, in those 
scenarios too. 
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APPENDIX: THE POWER SPECTRUM OF 5tf 

The entropy perturbation is given by the square of the 
field perturbation, and so we need the spectrum of that 
quantity. We assume that the power spectrum of 8<p itself 
is given by 



V Sl p(k) = A 



^H c a c / 

The spectrum of S(f> 2 is given by [p"8|f| 



(Al) 



k 3 r° H *^(|p|)^(|k- P |K 3 



2^ Jo 



| P | 3 |k-p|3 



d 6 p, (A2) 



where |p| and |k — p| are the wavenumbers. This inte- 
gral needs to be carried out with some care concerning 



the angular part. We can use rotational symmetry to 
align k along the z-axis to make |k — p| independent of 
the azimuthal angle, so that d 3 p = 2irp 2 sin 9 d6 dp. Sub- 
stituting in the spectrum from Eq. (Al) gives 



3 \2 



k A 



k -p|" sin 9d6 dp, 
(A3) 



where p = |p| and k = |k|. Using 



|k- P | 2 Hk| : 
leads to the integral 



2|k| |p|cos0, 



(A4) 



?V(fc) = 



k 2 A 2 



1 



(a e H e ) 2n 1 - n 
P n ~ 2 [\k-p 





Finally, setting u = p/k gives 

i- ' /■ \ 2n 
T^s^ 2 {k) = —— — ( — -TT- \ - 



1 — ri \a H c 

a c H c /k 



(A5) 

- {k+p) n - 1 } dp. 

(A6) 



[II 



1 



du . 



where we are interested in n positive, and usually less 
than one. This final integral has interesting properties. 
It is convergent as u ~ 0. For n < 1 the integrand is di- 
vergent at u = 1 (i.e. p — k) but it is integrable. Finally, 
at large u the integrand goes as u 2n ~ 4 . If n > 3/2, the 
upper limit dominates and one finds Vg^ cx k 3 ; this limit 
appears for instance in the analysis of preheating mod- 
els in Ref. |l9| where n = 3. However for the smaller n 
values we are interested in, the integral is dominated by 
U *~ 1. The upper limit of the integral becomes irrelevant 
and can be taken to infinity for modes with k < a c H , 
so that the integral is independent of k. The spectrum 
therefore has tilt Vs<j> 2 (k) <x k 2n . For the range of n we 
will be interested in, < 2n < 1, we can write 



V S4 ?{k) 



A 2 I{n) 



a c H r 



(A7) 



where 



I(n) = 



II 



iY 1 - 1 - (1 + uY 1 - 1 ] du. (A8) 



We are not aware of an analytical evaluation of this in- 
tegral except for special cases, but it is readily done nu- 
merically and we show the result in Fig. 1. 



tin Ref. |18) an unusual mix of Fourier conventions is used, 
but the final result is correct. 
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